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Abstract 

We present electronic and magnetic properties of a honeycomb compound In3Cu2V09 in this 
paper. We find that undoped parent phase is a charge transfer insulator with an energy gap of 1.6 
eV. Singly occupied 3z^-r^ electrons of coppers, similar to the pz electrons of carbons in graphene, 
lie from -5 to -7 eV and contributes an antiferromagnetic moment with S=l/2. Nonmagnetic 
vanadium, out of the expectation, has two electrons in V 3d orbitals. Oxygen 2p orbitals hybridizing 
with a small fraction of Cu 3d orbitals dominate the density of states near Ej?. We also estimate that 
the planar nearest-neighbor and next-nearest-neighbor superexchange couplings of Cu spins are 
Jiwl7.5 meV and J2~0.76 meV, giving rise to a low-dimensional antiferromagnet ^. We suggest 
that hole doping on O sites may bring about carriers and considerably enhance the screening, 
driving the system from an antiferromagnetic state to a spin liquid, or a superconducting ones, 
very similar to cuprates. 
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I. INTRODUCTION 



Electrons in reduced dimensional lattice, for example, in two-dimensional honeycomb 
lattice, has less coordinate number and small spin number, they may experience strong spin 
fluctuations and probably form extremely low-dimensional antiferromagnet, or even exotic 
spin liquid state, or quantum spin Hall state, so it has attracted great interests in recent 



n 



theories [ll? -[l2|. A recent synthesized honeycomb compound In3Cu2V09 [l4j-ll6| is such a 
system with reduced dimensionality. According to the chemical valence analysis, the copper 
ions with 3d^ conflguration contribute a spin-1/2 magnetic moment, while the 3d orbitals in 
vanadium ions are empty and contribute no magnetic moment. Due to the large separation 
between Cu/V-0 layers, the copper spins form a quasi-two-dimensional honeycomb lattice 



15| . This compound exhibits unusual magnetic behavior: it shows neither an AFM long- 



range order in magnetic susceptibility measurement, nor any magnetic phase transition peak 

in speciflc heat data [l^ [igI . This arises great interesting whether the ground state of this 

compound is an example of low- dimensional antiferromagnt (AFM), or a concrete example 

n 

of spin liquid system predicted by Meng et al. [1]. 

Additionally, it is well known that the electrons in graphene is a typical example of 
weakly interacting relativistic quantum systems. It is not very clear that how the properties 
of Dirac particles evolve with the increase of Coulomb interactions. Therefore the electronic 
motions in present transition- metal oxides In3Cu2V09 with the honeycomb structure is 
possibly a realistic testifying fleld of relativistic quantum mechanism for testifying how the 
electronic correlations or strong interaction affect the the motions and properties of the 
relativistic Dirac particles. In the mean time, the magnetic susceptibility and speciflc heat 
measurements are hard to distinguish possible ground state from a low- dimensional AFM 
and a spin liquid state 



16]. All of these drive us to explore the electronic properties and 
magnetism in In3Cu2V09 by combining the flrst-principles electronic structures calculations 
and analytic methods, since it could provide us some flrst insight to this compound. 

Utilizing the local density functional (LDA) and the correlation correction approaches, 
together with the analytical perturbation method, we show that undoped In3Cu2V09 is 
an AFM charge transfer insulator with copper 3z^-r^ orbitals half fllled and contributing 
S=l/2 spins. Out of our expectation, 3d orbitals in nonmagnetic vanadium have about two 
electrons. The oxygen 2p orbitals dominate the density of states (DOS) near E^. We not 



only extract tight-binding parameters between copper ions, but also show that In3Cu2V09 
falls in the two-dimensional AFM regime. We also suggest the possibility that doping drives 
an AFM-spin liquid or superconductor transition. The rest of this paper is organized as 
follows: in Sec. II we briefly describe the unit cell of In3Cu2V09 and numerical calculation 
methods; in Sec. Ill we present the major numerical results and our analysis; Sec. IV is 
devoted to the conclusion. 



II. NUMERICAL METHODS 

We first study the electronic properties of the honeycomb compound In3Cu2V09 by em- 
ploying the first-principles electronic structure calculation approach so as to elucidate the 
groundstate properties of In3Cu2V09. With the structural data measured by neutron pow- 



der diffraction we selected 30 atoms as a magnetic unit cell In6Cu4V20i8, as seen in 
Fig.l. The electronic structure calculations were performed using the self-consistent full po- 
tential linearized augmented-plane-wave (FP-LAPW) scheme in the WIEN2K programme 



package. [17|. We used 36 k points in the irreducible part of the first Brillouin zone. The 
muffin-tin sphere radii were selected to be 2.36, 1.93, and 1.71 a.u. for In, Cu (V), and O in 
In3Cu2V09. The plane-wave cutoff parameter RMT^max was 7.0, and the cutoff between the 
core and valence states was -6.0 Ry in all calculations. Exchange and correlation effects were 
taken into account in the generalized gradient approximation (GGA) by Perdew, Burk, and 



Ernzerhof (PEE) 18|. In order to take into account more explicitly the correlated motion 
of the 3d electrons of Cu (V) ions, we also performed GGA+U calculations for In3Cu2V09. 
The effective Coulomb interaction Ueff = U-J/^- (U and Jh are the on-site Coulomb inter- 
action and the Hund's rule coupling, respectively) was used instead of U. We take Ue//=8.0 
eV for copper and 4.0 eV for vanadium. 

To obtain the tight-binding parameters to model the low energy process in In3Cu2VOg 
in the paramagnetic phase, we also perform the local density functional (LDA) calcula- 
tions. The paramagnetic band structure calculation of In3Cu2V09 using the full-potential 
linearized augmented plane-wave plus local orbitals (FP-LAPW+lo) scheme implemented in 



the WIEN2K package IT]]. In order to compare our numerical results with the experimental 
data, we adopt the experimental structural data of In3Cu2V09 [14]. The parent material 
In3Cu2V09 has a honeycomb layered structure with Cu atoms forming hexagonal net and 




FIG. 1. (Color online) (a) Magnetic unit cell of In3Cu2V09, containing two crystallography primary 



V atoms in the center of an hexagon. The low-temperature experimental lattice parameters 



are a=3.3509A and c=11.9012A 



14l |. In each unit cell there are four Cu atoms, and each Cu 



is coordinated by five O atoms. Detail paramagnetic band structure results are addressed 
in what follows. 



A. Paramagnetic State 



We first present the paramagnetic {i.e. nonmagnetic) electronic structures of In3Cu2V09 
in Fig. 2, which is obtained in the FP-LAPW scheme through the WIEN2K package. It shows 
that the major orbital character of the band structures near Ej? is the — symmetry 
of copper 3d orbitals. This Cu 3d322_r2 orbital in the honeycomb lattice, similar to the p^ 
orbital of carbon in graphene, contributes an outstanding property of the band structures 
in Fig.2, i.e. around the H point a Dime cone with the approximately linear spectrum 
is observed. The linear energy spectrum ranges from -0.2 to 0.2 eV. We notice that the 
original Dime point around the K points opens a small energy gap, about 0.2 eV, also as 
seen in Fig.2, which is attributed to the weak interlayer coupling between Cu/V-0 hexagonal 
planes. Though the part of empty V 3d orbital are near the top of 3d3^2_^2 orbitals, they 
do not hybridize with each other. A part of filled oxygen 2p orbitals are close to the bottom 
of 3d322_r2 orbitals, these 2p orbitals also do not mix with 3d322_r2 orbitals, as we see in 
Fig.2. Thus the low-energy physics in In3Cu2V09 could be approximately described by the 
Cu 3d3^2_r2 orbitals across Yip. 

According to the Cu 3d322_r2 orbitals mainly distributed from -1.0 eV to 0.4 eV, we 
extract these four Cu 3d322_^2 bands from the full band structures and fit these four bands 
with a single-orbital tight-bind model with the nearest-neighbor (NN), the next NN and the 
3r(i NN hopping integrals in the Cu/V-0 plane and between the planes. These tight-binding 
parameters are listed in Table 1. One can easily find that the intraplane NN hopping integral 
ti is about 6 and 4 times of the intraplane NNN hopping and the interplane NN hopping, 
respectively, and dominates the transfer process in In3Cu2V09. 

5 




FIG. 2. (Color online) Band structures of In3Cu2V09 in the paramagnetic situation in the LDA 
framework. 
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B. Antiferromagnetic Ground State 

To understand the insulating nature of undoped compound, we further adopt the 
correlation-corrected GGA (GGA+U) scheme to uncover the groundstate properties of 
strongly correlated ln3Cu2V09. We consider four kinds of different magnetic configurations, 
nonmagnet, ferromagnet (FM), NeeJ AFM and A-type AFM, i.e., interlayer AFM and in- 
tralayer FM), as the possible initial states. Our numerical results show that the total energy 
of the nonmagnetic state is the highest, and that of the Neel AFM state is the lowest. The 
FM and A-type AFM configurations are almost degenerate. The energy difference between 
the lowest Neel AFM and the second lowest A-type AFM or FM states is about 58.5 meV, 
suggesting that the magnetic coupling between Cu spins in ln3Cu2V09 is dominantly AFM. 

The system is obviously an insulator with an energy gap of 1.6 eV, as seen in Fig. 3. 
After projecting various bands into each atom, we find that oxygen 2p orbitals consist of 
the major bands near E^?, ranging from to -4.0 eV. While the gravity center of Cu 3d 
orbitals lie far from Yjp. As shown in Fig. 3, the center of the spin-up 3d levels hes at 
about -5.9 cV, and that of the spin-down 3d orbitals lies at about -4.8 eV. This gives rise 
to the exchange splitting of Cu 3d orbitals about Ae^^l.! eV, implying that the Hund's 
rule couphng ^H—^exf^S^ 1.1 eV. Such a large exchange spfitting also shows the strong 
Coulomb correlation on Cu 3d orbitals. Notice that the charge transfer energy between O 
2p orbital and Cu 3d orbital in In3Cu2V09, A.—e2p — e^a, which is the same order to the 
insulating energy gap, is considerably smaller than the on-site Coulomb correlation of Cu 
3d electrons, U ~ Ue//= 6.0-8.0 eV, demonstrating that undoped In3Cu2V09 is a charge 
transfer insulator. 

TABLE I. Tight-binding parameters and superexchange spin coupling strengths Jq, in 
In3Cu2V09. Subscripts 1,2,3 and c represent the NN, next NN, 3rd NN hopping integrals or 
spin coupling strengths in the planes, and the NN ones between the planes, respectively. All the 
energies are measured in units of meV. 

tl t2 t3 tc Ji J2 Jc 

-180.2 -24.1 5.8 -48.9 16.2 0.3 1.2 
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FIG. 3. (Color online) Band structures of In3Cu2V09 within the spin-polarized GGA+U approach. 
The bands marked with heavy circles are Cu 3d orbitals for spin down channel. 
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The charge distribution and occupation number in the 3d orbitals of Cu and V are most 
essential. Our results show that in the Cu 3d orbitals, the hole mainly occupies the active 
d322_j,2 orbital. It has a total number of 1.17 electrons, or 0.83 hole, in the d322_^2 orbital. 
However, out of the early expectation and in contrary to the chemical valence analysis, V has 
totally 1.87 electrons in its 3d orbitals, though it shows nonmagnetic or no spin polarization 
in the present GGA+U ansatz. This arises from a strong p — d hybridization between V and 
O atoms, a fraction of electrons transfer from oxygen 2p orbitals to vanadium 3d orbitals. 



C. Density of States (DOS) 

The density of states (DOS) projected onto each atom in Fig.4 clearly demonstrates the 
level distributions of Cu, V, O and In atoms. Fig. 4a displays that the main DOS of copper 
lies from -7 eV to -4eV, though a tiny portion of 3d electrons lie between -4 eV to E^? which 
arises from the hybridization of Cu 3d orbital with the 2p orbitals near of surrounding 
O or with the 3d orbitals of surrounding V atoms. One sees from Fig. 4a that the 3d322_r2 
orbital in Cu is almost fully spin polarized. The magnetic moment of Cu spin is contributed 
from the well localized spin-down ?)d^z'^_^2 orbital at the position of about 7 eV below 'Ep. 
The sharp DOS peak of Cu in Fig. 4a is consistent with the spin-down flat band of Cu in 
Fig.3. Ourorbital ground state is the d3^2^^2 character, in agreement with the electron spin 



resonance experiment by Kataev et al. 19|. 



Meanwhile, small but finite DOS distributed from -6 eV to E^? is seen in V 3d orbitals, 
showing that the two 3d electrons form a wide energy band below the Fermi energy. This 
implies that it not only hybrids with the nearest-neighbor oxygen, but also mixes with 
distant Cu 3d electrons. Obviously V 3d orbitals contribute the lowest empty conduction 
bands in Fig.3, as we can see in Fig.4b. The projection of the DOS to each 3d orbital of V 
shows that the two 3d electrons in V are not spin polarized and allocates in five 3d orbitals 
with almost equal weight, in sharply contrast to those of Cu atoms. Fig. 4c shows that In 4p 
orbitals give rise to a small DOS in a wide energy range. In a contrast, the DOS of oxygen 
is greatly larger than that of In and mainly distributes near 'Ep, as seen in Fig.4d. Such a 
strongly correlated character in In3Cu2V09 is very similar to that in the parent phases of 
high-Tc cuprates. This implies that the doping in O sites can readily affect the transport 
properties in In3Cu2V09, and easily to drive the system transit to superconducting phase 
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FIG. 4. (Color online) Partial density of states (DOS) of In3Cu2V09 within the spin-polarized 
GGA+U method, (a). Cu 3d orbitals; (b) V 3d orbitals; (c) In 4p orbitals, and (d)0 2p orbital at 
site 13. 

once hole carriers are doped in oxygen 2p bands. This also suggests that doped In3Cu2V09 
will resemble to doped cuprates in many respects. 



D. Magnetism and Magnetic Couplings 



The groundstate magnetic properties of In3Cu2V09 is the central issue of this paper. As 
described in the above, our numerical results demonstrate that the stablest phase is Neel 
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AFM. The magnetic moment mainly contributes from copper spins with a net total magnetic 
moment of 0.78 fiB per Cu atom, in which the average electron number is 0.97 in the spin- 
up band, and 0.19 in the spin-down one. The magnetic moment mainly arises from the 
d^z^-r^ orbital of copper. Apparently, a strong electron correlation and a large Hund's rule 
coupling contribute such a large spin polarization, leading to the well defined local moment 
in coppers. In contrast, vanadium is nonmagnetic, though it is 3d^ configuration with two 
electrons distributing from E^? to -7.0 eV. Its nonmagnetic character may arise from the 
strong pd hybridization between V and O, resulting in a weak correlation in V 3d orbitals. 
We do not rule out the possibility of spin frustrations in the present Cu-V geometry, this 
may deserve further discussion in the future. Due to the full fillings of the oxygen 2p orbitals 
and In 4p orbitals, both O and In are not spin polarized. 

The insulating nature in the undoped In3Cu2V09 suggests that the local spins of coppers 
interact through a superexchange couplings mediated via oxygens and vanadium. One ex- 
pects that the possible contributions to the NN and next NN Cu-Cu superexchange couplings 
arise from the direct hopping between d^z^^r^ electrons, the indirect hopping between Cu 
spins through the pd hybridization of intermediate oxygen anions, and the indirect hopping 
through the 0-V-O bridge. To quantitatively obtain the superxchange coupling strengths, 
we utilize the tight-binding parameters obtained in Table I to estimate the NN and next NN 
superexchange couplings between Cu spins in the Cu/V-0 plane and between the planes, 
taking the effective Coulomb interaction U as 6 eV. We obtain that the intraplane and inter- 
plane NN superexchange couplings are about Ji ^16.2 meV and Jc ~1.2 meV, respectively, 
and the intraplane next NN superexchange couplings are about J2 ~0.3 meV. More long- 
range spin couplings are so small to be negligible. This gives rise to J2/J1 ~ 0.02, falling in 
the parameter range of a two-dimensional Neel AFM in the J1-J2 Heisenberg model j^, I22I . 
in addition to a rather weak interlayer coupling of Jc/Ji~0.07, suggesting strong quasi- 2- 
dimensionality. All of these values have been listed in Table I. Therefore, we conclude the 
undoped In3Cu2V09 is a strongly two-dimensional AFM. 

We notice that the critical value between the low-dimensional AFM phase and the spin 

n 

liquid one in the Hubbard model [51 is considerably smaller than that of the J2-J1 Heisenberg 



model |22l], we attribute to that the coexistence of charge and spin fluctuations expands the 
spin liquid or disorder region in the Hubbard model, in comparison with the Heisenberg 
model with only the spin fluctuation. From our analysis to the electronic state undoped 
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In3Cu2V09 can be well described by the J2-J1 Heisenberg model, with the Ji, J2 and Jc 
listed in Table I. 



III. REMARKS AND CONCLUSIONS 

Since the experimental data of both the specific heat and magnetic susceptibility do not 
exhibit AFM order-paramagnet transition 16||, the specific heat is almost dependent, 



and the magnetic susceptibility linearly increases with the lift of temperature, one natu- 
rally suspects that the ground state of undoped In3Cu2V09 might be a gapless spin liquid 
phase, a long-time searched exotic quantum phase both theoretically and experimentally, 
since a gapless spin liquid phase in Kagome lattice can also demonstrate similar unusual 



temperature-dependent behavior in the low T 



21 



20|. In the mean time a strongly low- 
Our theoretical analysis based on the 



dimensional AFM may exhibit similar behaviors 
first-principles electronic structures calculations and band structures fitting support the lat- 
ter, though it is not easily to understand why the system does not show Curie- Weiss law in 



the high temperature regime IGj. This exotic property is distinctly different from undoped 
AFM cuprates. 

Moreover, as we point out above, the hole doping in O sites of In3Cu2V09 may introduce 
carriers near Fermi energy. With the increase of carriers, one expects that the screening 
effect can greatly reduce the Coulomb interaction, and weaken the electronic correlation. 
This may drive the In3Cu2V09 system from a Neel AFM ground state to a spin liquid 
phase jl!l, since in the intermediate correlation regime a half-filled Hubbard system with the 

nn 

honeycomb structure exist a spin-liquid ground state [1|, Nevertheless, we also expect 
that doped In3Cu2V09 will not form the well-known Zhang- Rice singlet in the copper-oxygen 
plane of doped cuprates, since the latter arises from the hybridization between Cu ?)(\x'^-y2 
and O 2p orbitals, forming delocalized states in the xy plane; while in the former the 3d322_r2 
orbital mixes with O 2p orbital, forming delocalized state in the z-direction. This results 
in completely different transport properties in doped compounds, as Yan et al. observed in 



In3Cu2-xCoa;V09 and In3Cu2_xZn^.V09 [l6|. 



In summary, the electronic structure and properties investigation on quasi-two-dimensional 
honeycomb compound In3Cu2V09 demonstrate that undoped phase is a charge transfer in- 
sulator with a gap of 1.6 eV , and oxygen 2p orbitals dominate the bands near E^?. Spin- 1/2 
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local moment antiferromagnetically interacts and occupies the 3z^-r^ orbital in copper. In 
contrast, two electrons lie in the 3d orbitals of V and are nonmagnetic. The tight-binding 
parameters and estimated spin coupling strengthes suggest that undoped In3Cu2V09 is a 
strongly two-dimensional Neel antiferromagnet. Hole doping in oxygen sites may provide 
the testify field of the strongly interacting Dirac particles, suggesting an prospective future 
in doped In3Cu2V09. 



ACKNOWLEDGMENTS 

We acknowledge X. H. Chen provide us original experimental data. This work was sup- 
ported by the grants of the NSFC of China and of the Chinese Academy of Sciences. Nu- 
merical calculations were performed at the Center for Computational Science of CASHIPS. 



[1] Z. Y. Meng, T. C. Lang, S. Wessel, F. F. Assaad, and A. Muramatsu, Nature 464, 847 (2010). 

[2] S. M. Yan, D. A. Huse, and S. W. White, Science 332, 1173 (2011). 

[3] S. Saremi and P. A. Lee, Phys. Rev. B 75, 165110 (2007). 

[4] A. Vaezi and X.-G. Wen, larXiv:1010:57441 larXiv: 1101. 16621 

[5] B. K. Clark, D.A.Abanin, and S. L. S ondhi, Phys. Rev. Lett. 107, 087204 (2011). 

[6] F. Wang, Phys. Rev. B 82, 024419 (2010). 

[7] A. Liebsh, Phy Rev. B 83, 035113 (2011). 

[8] J. He, S.-P. Kou, Y. Liang and S.-P. Feng, Phys. Rev. B 83 205116 (2011); ihid, 84, 035127 
(2011). 

[9] Y.-M. Lu and Y. Ran, Phys. Rev. B 84, 024420 (2011). 

[10] H. H. Zhao, Q. N. Chen, Z. C. Cai, M. P. Qin, G. M. Zhang, and T. Xiang, arXiv: 1105. 27161 

[11] W.-S. Wang, Y.-Y. Xiang, Q.-H. Wang, F. Wang, F. Yang, and D.-H. Lee, Phys. Rev. B 85, 
035412 (2012). 

[12] W. Wu, S. Rachel, W.-M. Liu, and K. Le Hur. larXiv: 1106.09431 

[13] T. Li, Euro. Phys. Lett. 97, 37001 (2012). 

13 



[14] A. Moller, U. Low, T. Taetz, M. Kriener, G. Andre, F. Damay, O. Heyer, M. Braden, and J. 

A. Mydosh, Phys. Rev. B 78, 024420 (2008). 
[15] M. Yehia, E. Vavilova, A. Moeller, T. Taetz, U. Loew, R. Klingeler, V. Kataev, and B. 

Buechner, Phys. Rev. B 81, 060414(R) (2010). 
[16] Y.-J. Yan, Z.-Y. Li, T. Zhang, X.-G. Luo, G.-J. Ye, Z.-J. Xiang, P. Cheng, L.-J. Zou, and X. 

H. Chen. [arXiv:1106J7l3l accepted to Phys. Rev. B (2012). 
[17] P. Blaha, K. Schwarz, G. Madsen, D. Kvasnicka, and J. Luitz: "Computer Code WIEN2k, an 

augmented plane wave plus local orbitals program for calculating crystal properties", Karl- 

heinz Schwarz, Technische Universitst Wien, Austria, (2001). 
[18] J. P. Perdew, K. Burk, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996). 
[19] V. Kataev, A. Moller, U. Low, W. Jung, N. Schittner, M. Kriener, and A. Freimuth, J. Magn. 

Magn. Mater. 290-291, 310 (2005). 
[20] Y. Ran, M. Hermele, P. A. Lee, X.-G. Wen, Phys. Rev. Lett. 98, 117205 (2007). 
[21] K. Oka, L Yamada, M. Azuma, S. Takeshlta, K. H. Satoh, A. Koda, R. Kadono, M. Takano, 

and Y. Shlmakawa, Inorg. Chem, 47, 7355 (2008). 
[22] F. Mezzacapo and M. Boninsegni, Phys. Rev, 85, 060402 (2012). 



14 



